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Qh Abstract: We explore, in a concrete example, to which extent new particle mass deter- 

minations are practicable with LHC data. Our chosen example is that of Yukawa-unified 
SUSY GUTs, whose viability has been recently studied for two general patterns of soft 
SUSY- breaking terms. We note that both patterns of SUSY spectra do not admit long 
decay chains, which would make it possible to determine the masses of the SUSY particles 



involved using endpoints or mass relations. We thus take the so-called mT2-kink method 
as our key strategy, since it does not rely on the presence of long decay chains. We then 
discuss a procedure allowing to determine the masses of the gluino, of the lightest chargino 
as well as of the first two neutralinos and, for the scenario where a stop is lighter than 
the gluino, the mass of the light stop too. Our worked example of Yukawa-unified SUSY 



GUTs may offer a useful playground for dealing with other theories which predict similar 
patterns of SUSY spectra. 
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1. Introduction 

Existing data on established collider quantities - in particular electroweak precision tests, 
quark masses and flavour-changing neutral current processes - provide crucial constraints 
on many sets of Standard Model (SM) extensions, whose new interactions do in general 
imply tree- or loop-level deviations in some or all of these observables. These tests allow 
to learn whether a given class of theories is viable at all, and, if so, to learn about the 
features of the viable parameter space. However these tests, at present, are mostly null 
tests, where the chosen class of theories is required to produce a small signal as compared 
with the SM. It is clear that, in order to single out the considered class of theories, one has 
to find observables for which an unequivocally different behavior with respect to the SM is 
predicted. With the imminent flow of direct data from the Large Hadron Collider (LHC), 
such observables may be found in inclusive searches or, better, in the direct determination 
of at least the lightest part of the new particles' spectrum. The second possibility is clearly 
more ambitious, but also much more powerful than the first one to tell apart a given theory 
against other possibilities. Aim of the present paper is to explore, in a concrete example, 
to which extent new particles' mass determinations are practicable with LHC data. 

The example we focus on is that of supersymmetric (SUSY) grand unified theories 
(GUTs) with Yukawa unification (YU), whose viability and expected signatures have been 



extensively studied with different approaches for various general patterns of soft SUSY- 
breaking terms [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19]. x We will focus 
here on the two scenarios considered in refs. [10] and [16], that we briefly summarize in 
the rest of this section. In ref. [10] the case where soft-breaking terms for sfermions and 
gauginos are universal at the GUT scale was considered. It was found 

(a) that this class of theories is phenomenologically viable only by advocating partial 
decoupling of the sfermion sector, the lightest mass exceeding 1 TeV; 

(b) that phenomenological viability can be recovered without decoupling by relaxing 
t — b — T unification to b — r unification, equivalent to a parametric departure of tan j3 
from the value implied by exact YU. This solution is non-trivial since the constraints 
from mi, and respectively FCNCs prefer high, O(50), and respectively low values of 
tan/3. Indeed, a compromise solution has been found to exist only for the narrow 
range 46 < tan/3 < 48, implying that the breaking of t — b YU must be moderate, in 
the range 10 - 20%; 

(c) that the requirement of b — r unification and the FCNC constraints are enough to 
make the predictions for the lightest part of the SUSY spectrum robust ones in the 
interesting region of point (b). These include a lightest stop mass > 800 GeV, a gluino 
mass of about 400 GeV and lightest Higgs, neutralino and chargino close to the lower 
bounds. Specifically, the lightest neutralino is in the ballpark of 60 GeV, and the 
lightest chargino about twice as heavy. The rest of the SUSY spectrum lies instead 
in the multi-TeV regime: in particular first and second generation sfermion masses, 
that almost do not RGE-evolve from the GUT scale to the Fermi scale, remain of 
the order of the GUT-scale bilinear mass m\s > 7 TeV. For an illustrative example 
of this spectrum, see the leftmost column of table 1 . 

Ref. [16], instead, explored the complementary case where YU is kept exact, while 
the requirement of universal soft terms at the GUT scale (not very justified from a theo- 
retical standpoint) is relaxed. Ref. [16] focused on the scenario in which the breaking of 
universality inherits from the Yukawa couplings, i.e. is of minimally flavor violating type 
[32, 33]. This general setup allows in particular for a splitting between the up- type and 
the down-type soft trilinear couplings. Comparison with data revealed for this trilinear 
splitting scenario the following features [16]: 

(i) agreement with data clearly selects the parametric region with a large /i-term and 
a sizable splitting between the ^4-terms. At the price of a slight increase in the fine 
tuning required to achieve EWSB with precisely the correct value of Mz, this scenario 
allows a substantial improvement on other observables that, on a model-dependent 
basis, do often require some amount of fine tuning as well; 

(ii) in particular, and quite remarkably, phenomenological viability does not invoke a 
partial decoupling of the sparticle spectrum, as is the case in the scenario of ref. 



^or earlier literature on the topic of YU within SUSY GUTs, the reader is referred to [20, 21, 22, 23, 
24, 25, 26, 27, 28, 29, 30, 31]. 
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Table 1: Spectrum predictions for representative fits of the viable scenarios studied in refs. [10] 
and [16]. All masses are in units of GeV. Uppercase and lowercase masses stand for pole and 
respectively DR masses. 



[10], but instead it requires various SUSY particles, notably the lightest stop and the 
gluino, to be very close to their current experimental limits. Specifically, the lightest 
stop is lighter than the gluino, in the range [100, 200] GeV. This testable difference 
reflects the substantial difference in the mechanism that makes the two classes of 
models phenomenologically viable. Ultimately, it maps onto the difference in the 
assumed form for the GUT-scale soft terms; 

(Hi) concerning the rest of the spectrum, the predictions for the gluino, the lightest 
chargino and neutralino masses, as well as those for the heavy part of the SUSY 
spectrum, are very similar to the corresponding ones in the scenario of [10] (see ref. 
[16] for further details). Again, because of the cross-fire between the many constraints 
used, all these spectrum predictions, including those of point (ii), are robust ones. 



As already mentioned, this paper is devoted to exploring the possibility of directly 
measuring the masses of the lightest SUSY particles in the scenarios of refs. [10, 16]. In 
particular, it will discuss a strategy - applicable already to data collected at the LHC - that 
permits the determination of the masses of the gluino, the lightest chargino and neutralino 
and, for the scenario of ref. [16], the lightest stop also. This task clearly requires to focus 
on short decay chains, and makes extensive use of so-called "transverse mass" variables. 

Our strategy will be described in section 3, after an introduction to the methods for 
new particles' mass determination at hadron colliders, to be presented in the next section. 



2. Methods for SUSY mass determinations at the LHC 

In this section we briefly recall the methods devised so far for new-particles' mass deter- 
minations at the LHC, and discuss their relation with our problem of interest. A compre- 
hensive and very recent review can be found in ref. [34]. 

2.1 Methods other than mj2 

The most straightforward mass-determination method is, in general, to find a peak in the 
invariant mass distribution of the decay products of the new particle of interest. This is 
however inapplicable if the final state includes a component that escapes detection, which 
is our case, because of the missing 'lightest SUSY particle' (LSP). 

For decay chains that include LSPs, in order to reconstruct SUSY masses one can 
exploit kinematic relations involving these masses and the (measured) momenta of the vis- 
ible particles. Among the strategies devised to this end, a first possibility is given by the 
so-called endpoint method [35, 36, 37, 38, 39, 40, 41, 42, 43]. It is based on the observation 
that, in a given decay chain, the endpoint values of the invariant mass distributions con- 
structed for visible decay products depend on the masses of the invisible particles as well. 
A prototype example is provided by the squark decay chain [44] 

Q -»• X2Q ->• PFq ->• Xit + rq , (2.1) 

which would be available for nig > m^o > m» > m^o . Assuming that the four-momenta of 
the leptons and of the g-initiated jet can all be determined with reasonable accuracy, the 
above decay process provides four invariant-mass distributions, i.e. those of nia, TW(high) = 
max(m q £+,m q £-), m q £(\ ow ) = min(m^+, m q £-), and m q u- Then their endpoints can be 
inverted to give all the four sparticle masses involved in the decay. However this decay 
mode is kinematically closed for our scenarios of interest, since all the sleptons are (much) 
heavier than any of the other sparticles involved in eq. (2.1). 

Besides, it can be shown that, in order for all the sparticle masses to be, even in 
principle, separately reconstructible, the endpoint method requires necessarily a long decay 
chain with at least 3 decay steps [45]. If this condition is not fulfilled, the number of 
measurable invariant masses is smaller than the number of unknown particle masses, and 
one can at best determine some combinations of mass differences, rather than determining 
the absolute mass values. A simple example would be a single-step 3-body decay 

X° 2 -+ x\t+t~ . (2.2) 

In this case, the dilepton invariant mass has an endpoint equal to m^o — m^o. In order 
to determine m^o and m^o separately, one would need additional kinematic variables as 
discussed in [44]. 

A second method is that of mass relations 2 [46, 47, 48, 49, 50]. The idea here is 
that the four-momenta of missing particles in the decay products can be reconstructed 
by exploiting various constraints including the on-shell constraints on the decay chain. 



2 This method is often called the polynomial method. 



Provided the number of branches in the decay chain is sufficiently high, the number of 
constraints can be made to exceed the number of unknowns, and one can solve (or actually 
fit) for all the unknowns (for a detailed discussion see [45]). For example, in an n-step 
cascade decay initiated by a SUSY particle Y, with each vertex branching into a SUSY 
particle Ij and a SM visible particle with (reconstructible) momentum pi , 

Y -► h + V(Pi) ~^h + V(p 2 ) + V( Pl ) -> ... -> I n (k) + V{p n ) + ... + U(pi) , (2.3) 

one has n + 1 independent mass-shell constraints, namely 

/c 2 =m\ ,(k+ Pn ) 2 = mj nl ,... ,(/c+pi + ...+p n ) 2 = m 2 , , (2.4) 

where k is the four- momentum of the missing particle I n . For N such events, the number 
of constraints will be N(n + 1), while the number of unknowns will instead be 47V" + (n + 1), 
i.e. the four-momentum of the missing particle in each event and the n + 1 masses of the 
SUSY particles which are common to all events. In order for the number of constraints to 
exceed the number of unknowns, one needs n > 4 and also N > (n + l)/(n — 3). 

A case where the mass-relation method would be more appropriate is that of sym- 
metric n-step decays of pair-produced SUSY particles, Y + Y, for which one can use the 
missing transverse momentum constraint also [48, 50]. In this case, again for iV events, the 
unknowns include as before n + 1 SUSY particle masses which are common to all events, 
and the four- momenta of the two missing particles in each event, so the total number 
of unknowns is given by (n + 1) + 8N. As for the available constraints, for each event 
one has 2(n + 1) constraints from the mass-shell conditions and two constraints from the 
missing transverse momentum, so 2N(n + 1) + 2N constraints in total. We then find the 
number of constraints is equal to or bigger than the number of unknowns if n > 3 and 
N > (n + l)/2(n — 2). These observations lead to the conclusion that one needs at least a 
3-step or a 4-step cascade decay in order for the mass-relation method to be applicable. 

On the other hand, for the scenarios of refs. [10, 16] the predicted SUSY spectra imply 
that, at the energies available at the LHC, there cannot be any long decay chain (n > 3) 
on which the above discussed endpoint or mass-relation methods might realistically be 
applied. Specifically, the only 3-step cascade decay with sizable number of events is the 
gluino decay in the scenario of [16] , g — >• tt\ — > bWbx± — > bqqbqqxo, which however suffers 
from very large jet combinatorics. As a result, the mass-relation method simply cannot be 
used in our case, and the endpoint method can determine at best mass differences. 

2.2 The mT2-kink method 

A third method exists, however, that is able to determine SUSY masses even in the absence 
of long decay chains. This method, called the rriT2-kink method [51, 52, 53, 54, 55, 56, 57], 3 
exploits the fact that the endpoint value of the transverse mass variable m,T2 > regarded as 
a function of the trial mass m x of the missing particle in the decay products, exhibits a 
kink at m x = m* ruc [51]. As the endpoint value of rriT2 at m x = m* ruo corresponds to the 



3 For another method, also applicable to short decay chains, see [58]. 



mother particle mass, the mT2-kink method determines both the mother particle mass and 
the missing particle mass simultaneously. As we will see, applying the mj^-kink method to 
pair-produced gluinos in the scenario of [16], while regarding charginos as missing particles, 
one can determine both the gluino mass and the chargino mass. Likewise, for the scenario 
of [10], one can consider the m^-kink of gluino pairs while regarding the second lightest 
neutralinos as missing particles. With those masses determined by the m^-kink, the rest of 
the light part of the SUSY spectrum can be determined by applying the endpoint method. 
The collider observable rriT2 is a generalization of the transverse mass rriT which has 
been introduced for decay processes producing invisible particles [59, 60]. Specifically, 
considering the decay Y — > V{p) + x(^)i where x(fc) is an invisible particle with four- 
momentum k, and V{p) stands for an arbitrary number of visible particles with total 
four-momentum p, the transverse mass my is defined as 



m T 



my + m 2 + 2 



m\r + \p T 



n m |2 



m 2 , + \k T \ 



p T ■ k T ) 



(2.5) 



where my and m x are the invariant masses of V and x, respectively. As the four-momentum 
of x cannot be directly measured, let k and m 2 , = k 2 denote a trial four-momentum and 
respectively a trial mass squared for x, which are introduced for the sake of discussion. A 
key feature of ttit is that it is bounded by the invariant mass, i.e. the true mother particle 
mass m Y nc , when m x and kx take the true values: 



m T (my 



true 



my 



+ (m 



m 



true \ 2 



kT 



£r me ) < (mTl 2 



+ 2 



m\ 



+ Wt\ 



( m true)2 + |fctrue|2 cos h( w 



^ Ue ) 



p T ■ 4 mG 



is the pseudo-rapidity 




Figure 1: Sketch of event topology rele- 
vant for the applicability of the my 2 -kink 
method. 



where rj = ^lnff^ 
While the true momentum of x cannot be di- 
rectly measured, its transverse component k T :ne 
can be inferred from the missing transverse mo- 
mentum if x is the only missing particle in the 
whole event: k^ ue = ft?. Then, once m* rue 
is known by some other information, one can 
determine my e from the endpoint value of mr(m x 
m£ uc , fc r = jfc£ uc =f T )- 

As typical SUSY events involve a pair of 
SUSY particles decaying into two invisible LSPs, 
the transverse mass should be accordingly gen- 
eralized. At the LHC, generic SUSY events take 
the form (see fig. 1) 

p + p^? U + Y + Y (2.6) 

->U + Vi(p)+x(k) + V 2 (q)+x(l), 



where U, often called the upstream momentum, denotes the total four-momentum of visible 
particles not coming from the decays of the SUSY particle pair Y + Y, x is the invisible 



LSP, and again each Vi (i = 1, 2) denotes an arbitrary number of visible particles produced 
by the decay of Y or of Y. (The upstream momentum might come from initial or final 
state radiation or from the decay of a heavier SUSY particle.) Unlike the case of single 
missing particle, the transverse momentum of each LSP cannot be determined in this case, 
although their sum can be read off from the missing transverse momentum. One can then 
introduce trial LSP transverse momenta, kr and It, under the constraint kx + W = Pt, 
and consider the transverse mass of each decay mode. Then the rriT2 variable for a trial 
LSP mass m x can be constructed as [59, 60] 

m T2 (event; m x ) = ming^^ [max (m T (Y -> Vi + x), m T (Y ->■ V 2 + x))] , (2.7) 

where each event is specified by the set of measured kinematic variables {my 1 , pr , fny 2 ,Qt,Pt}, 
with pr the total missing transverse momentum. Note that, due to the relation pr = 
~Pt — Qt — Ut, use of pr as an event variable is equivalent to using the upstream trans- 
verse momentum Ut- One can now observe two main properties of the vnri variable: (a) 
once the rest of the input is fixed (i.e. given the event), rriT2 is a monotonic (increasing) 
function of the trial LSP mass; (6) when the trial LSP mass is equal to the true LSP mass, 
m,T2 is bounded from above by the true mother-particle mass m'™. 



It follows that, studying mri as a function of m x in the neighborhood of m^ 
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See fig. 2 for an illustration. In fact, a relevant issue for the mT2-kink method 
is how sharp the kink is, which is equivalent to how much the slopes of mr2 curves (at 
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vary over the available event set. It has been noted that endpoint events with 
different my or different Ut generically have different slopes [51, 52, 53, 54, 55, 56]. 

This implies that the kink can be sharp enough 
to be visible if V corresponds to multi-particle 
states with wide range of my and/or if a large 
\Ut\ is available. As the possible range of my is 
bounded by m,y uc — m* me , the kink can be sharp 
enough only when m^ is substantially heav- 
ier than m* rue . As for \Ut\, although in prin- 
ciple it can be substantially larger than ?riy uc , 
in reality the best possible situation would be 
\U T \ ~m$ uc . 

2.3 Generalizations 

Although x was identified as the LSP in the 
above discussion, the very same argument ap- 
plies to the more general case where x lS a generic 
SUSY particle produced by the decay of Y or 
Y, as long as x(A;) and x(0 have the same mass 
[61, 45]. For instance, x might be a heavier neu- 




Figure 2: Pictorial representation of the 
rriT2 curve for two events (solid vs. bro- 
ken lines) in the endpoint region close to 
n r iT2(i r n x = m x Tno ) = my" . The x- and 
y-axes are in arbitrary units. 



tralino or a chargino, decaying into the LSP and some SM particles, provided one regards 



all those particles as missing particles. Indeed, as already mentioned, in the following we 
will be applying the m^-kink method to gluino-pair decay in the scenarios of [10, 16], with 
X being a chargino or the second lightest neutralino respectively 4 Looking at the SUSY 
spectra, summarized in table 1, the gluino undergoes a relatively short (2-step) cascade 
decay into x pl us several visible particles, and also the gluino mass is sensibly above m x . 
This by itself indicates that a visible kink might be possible even with events having a 
negligible \Ut\- 

In cases where the visible decay products involve multiple jets, it is often difficult 
to identify the endpoint position because of the smearing caused by poor jet momentum 
resolution. In such cases, to reduce the endpoint smearing, one can consider an alternative 
form of m,T2 defined solely in terms of the jet transverse momenta [51, 54]: 

m T2 = mm %T+rT= j T [max (m T (Y -»■ V x + x), m T (Y -> V 2 + x))] , (2.8) 

where now m/r(Y — > V + x) is given by 

ml = m 2 TV + m\ + 2 ( \jm\ Y + \p T \ 2 \jm\ + \k T \ 2 - p T • k T J . (2.9) 

Here 

m?t v (Y^V + x) = Y, m « + 2 Y,( E <* TE PT-PaT-P0T) , (2.10) 

a a >/3 

where p a T is the measured transverse momentum of the a-th jet in the set of visible decay 
products V, and E a T = \/ m a + |p q t| 2 , with m a denoting the a-th partem mass, rather 
than the jet invariant mass, i.e. m a = for non-6 jets and m a = rrib for 6-tagged jets. As the 
transverse momentum can be measured with better accuracy than the other components, 
the m,T2 defined in this way suffers less from jet momentum uncertainty. In addition, it 
still obeys all the features for a kink at m x = m* ruc with ?ri™| x (m x = m* mc ) = my at 
partonic level. On the other hand, as rriTV < m Vi the rriT2 defined with the transverse 
mass rriTV has generically smaller statistics near the endpoint than the mT2 defined with 
the full invariant mass my. Previous studies suggest that the reduced statistics is not a 
severe drawback compared to the gain from reduced smearing [51, 54], and thus the m,T2 
defined in terms of the jet transverse momenta can reveal a kink even when the kink of m,T2 
defined with my is smeared away due to poor jet momentum resolution. In the following 
analysis including detector effects, we will use the rriT2 defined with my if its distribution 
shows a clear endpoint, and the mx2 defined with m^y otherwise. 

To conclude this section, we would like to mention that several new methods of mass 
measurement for events involving invisible particles have recently been proposed, e.g. one 
exploiting the cusp structure of the distribution of certain kinematic variables [64] and 
another based on a generic algebraic singularity that arises in the observable phase space 
obtained by projecting out the unmeasurable kinematic variables [65]. It would be worth- 
while to investigate how much useful those methods may be for the specific SUSY scenarios 
of [10, 16] that we are studying in this paper. 

4 One can further generalize rriT2 to the case in which the parent or daughter particles are not identical 
[62, 63]. 
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Table 2: SUSY production cross sections from Pythia 6.42. Decay tables are calculated with 
SUSYHIT and then fed to Pythia. Missing entries are meant to be below 0.1 permil. 



3. Strategy 

In this section we apply the idea of mr2-kink described above to the SUSY GUT scenarios 
exemplified in table 1, and henceforth referred to as HS (heavy stop) scenario (ref. [10]) 
and LS (light stop) scenario (ref. [16]), because of the heavy and respectively light stop 
relatively to each other. A first necessary piece of information is that of the main 
SUSY production cross sections for the two scenarios. Table 2 shows the most important 
production mechanisms, for pp events with \/s equal to 10 or 14 TeV. 5 They have been 
calculated with Pythia 6.42 [66]. For both HS and LS scenarios the dominant SUSY- 
production mechanism is gg. Interestingly, for the LS scenario, tit* production is also 
large, close to 40%. Finally, also X2X1 associated production is non- negligible in both 
scenarios. The second needed piece of information is that of the main decay modes for 
the produced particles. These decay modes are reported in table 3, and calculated with 
SUSYHIT [67]. We will elaborate on these production and decay figures in due course 
during the analysis. 

In the light of the information in tables 2 and 3, the rest of this section is devoted to 
a short description of our mass-determination procedure, as it would be carried out in a 
parton-level analysis. In order to provide a quick overview of the various steps, the full 
procedure is also schematically reported in table 4. This table shows that our procedure 
is able to determine the g, Xi, X2> Xi masses in either scenario and, for the LS scenario, 
also the t\ mass. The concrete implementation of our strategy, carried out on 100/fb of 
LHC data at the design center-of-mass energy of 14 TeV, and including realistic detector 
effects, is postponed to section 4. 

3.1 Strategy: LS scenario 

Given the generally lighter masses and the higher cross sections for the LS scenario, our 



Our analysis will be focused on the 14 TeV case. The 10 TeV case is reported for guidance, with the 
aim of showing that production cross sections for all considered processes are such that our strategy should 
still apply at this value of y/s. 
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Table 3: Main decay modes in % for g, i\, x° an d Xi > calculated with SUSYHIT 



strategy is richer in this case, and we describe it first. 




Figure 3: Decay chain of interest for step 1 (LS 
scenario) . 



Step 1: Construct rriTi for gg —> 4 b + 

2 W (+ 2 £) + J)t decay events, thereby 
measuring m y and m~± 

In the scenario of ref. [16], the gluino de- 
cays always as g — > tit, with the top quark 
subsequently decaying as t — > bW and the 
t\ as Xi b- Requiring the Xi to decay as 
qq'xi would have the advantage that all the decay products apart from the x? would be 
reconstructible, in principle. However, in practice, the event trigger would consist of 2 b 
jets and 4 additional non-6 jets (2 jets from W and another 2 jets from \\ ) f° r each of 
the two decay chains, making this channel very challenging already on account of the jet 
combinatoric error. 

A more promising strategy is to focus instead on leptonically-decaying charginos, 
namely Xi ~~ *■ x\$- v li an d to regard the whole Xi -initiated decay chain (see fig. 3) as 
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Step 


Event trigger 


Event variable 


Allows to determine 



u 


1 


gg^Ab + 2W (+2£)+$ T 


m T 2 
(with 2£ g ^ T ) 


nig and m-± 

M Xi 


a 

CB 



2 


iiii -^2b + Aq + i/) T 


fnT,bqq 


l l Al 


IB 

1-J 


)) 


" 


1UT,qq 


Al Al 




3 


xtj&->t+e-? +i> T 


ma 


m^o — m^o 

A2 Al 



in 
IS 

d 

CB 



1 


gg -> 46 (+2 7 ) + ^ T 


m T 2 
(with 27 g $t) 


m§ and m^o 


2 


xfxg ->*+*-+*' + &. 


mu 


A2 Al 




3 


X^X2^2g + 2^ + ^ T 


m qq 


m~± — m^o 

Al A.1 




?? 


r> 


met 


A2 Al 



Table 4: Overview of our mass-determination strategy. 

•f>T- In this case the final state will contain 2b-\-W + £1/^ + x± f° r each decay chain, with £ 
either of e, [i. The ^t w iH then be the sum of the transverse momenta of the x\ an d of the 
Iv pair. Concerning the W, one can require it to decay in a qq' pair, so that m qq = m\y (up 
to the W width) in each of the two decay chains. The event topology to look at is therefore 
4b + 2W + 21 + $t'- while the combinatoric problem has not completely disappeared, it has 
been substantially mildened with respect to the case mentioned in the previous paragraph. 
Constructing the m-r^-kink for this event topology (with, as mentioned, the 2 leptons' 
Pt included in fa) allows to determine simultaneously rrig and fn~±. 



xi 



h 





Step 2: For t\t\ — > 2 b + 4 Q + fa events, deter- 
mine the mass differences mr — m,-o andm~±—m^a 

JJ *i Xi Xi Xi 

from the endpoints of the transverse-mass distribu- 
tions mxfiqg and mT, qq , respectively 

Figure 4: Decay chain of interest for In the scenario of ref. [16], given the lightness of the 
step 2 (LS scenario). t±, its pair-production is substantial and the only 

kinematically allowed decay mode is i\ — > X\ b- One 
can further require the X\ to decay (besides the invisible Xi) m to visible products (see 
fig. 4), which happens in about 20% of the cases. The mass differences m^ — m^o and 
m~± — m^o can then be obtained as the endpoints of the transverse-mass distributions 

Xi Xi ^ 

calculated respectively on the bqq' and qq' systems, for each of the two decay chains of the 
event. An evident problem here is to correctly assign the b- and ^-jets to the two decay 
chains. We will come back to this in the analysis. 

Step 3: For X\ X2 ~^ f- + ^~ + £' + fa events, determine the mass difference m^o — m^o 
from the endpoint of the invariant-mass distribution mu 
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Figure 5: Decay chain of interest for step 
3 (LS scenario). 



Within the LS scenario, information on the m^o 
mass can only be obtained from direct production 
channels, as a simple inspection of table 3 reveals. 
In particular, one can consider the quite large 
Xi X% associated production, and require both spar- 
ticles to decay to leptons (besides the Xi)-> namely 
xl -> x\^t~ and xf -> Xi^e (see fig. 5). Re- 
constructing the visible part of the x% decay allows 
to determine the mass difference m 



X2 



-m^o, hence 




Figure 6: Decay chain of interest for 
step 1 (HS scenario). 



rris-o, since the lightest neutralino mass is already known from the previous step 2. 

3.2 Strategy: HS scenario 

We next turn to discussing SUSY mass determinations within the HS scenario of table 1 . 

Step 1: Construct the mT2 variable for gg — > 4 b (+ 
2 "f) + j>T decay events, thereby measuring m g and 

777.^0 
X.2 

Gluino-gluino production, followed by gluino decay- 
ing as g —7- bl>X2, ^ s probably the golden mode for the 
HS scenario. If one further considers the radiative 
decay of the x% in to Xi (note that this decay mode 
has the largest branching ratio), the implied event trigger, namely 4b + 27 + p°T, has very 
small contamination from SM background. Including the gammas in the fa, the kink in 
the mT2 variable constructed for these events allows to simultaneously determine nig and 

m x° 2 - 

Alternatively, one might attempt to apply the m^-kink method to the subsystem 

involving 2x\ as parent particles and 27 as their visible decay products. In this case, 

the kink in mT2 arises due to nonzero upstream transverse momentum and it allows to 

determine the masses of x% and Xi- From the analysis we find however that the resulting 

kink is not as clean as the one obtained in the approach of the previous paragraph, where 

also the 46 are included as visible decay products. 



t' + i>T events, determine the mass difference m^o 



Step 2: For X f X2 -> e+ ?~ + 

from the endpoint of the invariant-mass distribution mu 

The implementation of this step is entirely analogous as in step 3 of the LS scenario. 



Xi 



in-o 



Step 3: For x.\ X2 ~ * %Q ~f~ %& ~t~ P*T events, determine the mass difference rn~± 

from the endpoint of the invariant-mass distribution m qq (and again m^o — m^o from mm) 

Still to be determined is the X\ mass. At variance with the LS scenario, this information 
cannot be easily extracted from secondary X\ generated in gg or i\t\ production events. 
However, one can consider X\ X2 associated production, with the X\ decaying as X\ ~^ 
XiQl' ■ Concerning the X21 ^ e radiative decay to Xi is, as already noted above, the one 
with the largest branching ratio. However, the resulting event topology, namely 2 q + 7 
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+ fop, is dominated by background events, that we were not able to get rid of. For this 
reason, we required the X2 to decay as y^U (and veto on £ = r, since its decay would 
produce missing momentum from neutrinos). The invariant-mass distribution of the qq' 
pair is bounded from above by m- 
from step 1. 






allowing to obtain m-±, since m^o is known 




In addition, the invariant-mass distribution 
met provides a cross-check with step 2 on the de- 
termination of m,-o 

A.2 



Xl 



4. Analysis and results 

Figure 7: Decay chain of interest for step l n this section we would like to show that the 
6 (Ha scenario). mass-determination strategy described in the pre- 

vious section is indeed practicable with real data 
collected at the Atlas [44] and CMS [68] experiments, and permits to determine most of the 
considered masses within about 20 GeV of the true values, with an integrated luminosity 
not exceeding 100 fb _1 at the design LHC center-of-mass energy of 14 TeV. To address this 
point, we have simulated full events with Pythia 6.42 [66] and modeled LHC detector effects 
with PGS4 [69] . Specifically, we adopted the default lhc . par card that comes with the 
latest PGS4 version. Purpose of this analysis is to get a presumably realistic idea of how 
much do the signals degrade once detector effects are taken into account. It goes without 
saying that the 'true' answer will require a full-fledged detector simulation, with insider 
knowledge of all acceptance and resolution details, which is prerogative of experimentalists 
only. 

We next proceed to the discussion of the analysis, following the steps described in the 
previous section. 

4.1 Analysis: LS scenario 

Step 1: Construct rriT2 for gg — > 4 b + 2 W (+ 2 £) + for decay events, thereby measuring 
ijig and m~± 

a Xi 

We select events with 4 jets tagged as b, 6 4 q jets (here and henceforth, we will so indicate 
jets not tagged as b) and 2 leptons. To reconstruct the W bosons, we make all the possible 
di-jet combinations out of the 4 q jets, and accept only events for which both di-jet invariant 
masses for at least one pairing are within My/ ± 20 GeV. In case that more than one paring 
fulfill this requirement, we choose the pairing with minimal 



/ AM^(2q)+AM^(2q'). 
Furthermore, in order to eliminate possible SM backgrounds, we impose the following cuts: 



(a) pn.2,3,4 > 50, 30, 20, 20 GeV for the 4 b jets; 

(b) Missing transverse energy $t> 50 GeV; 

For the sake of clarity, the fo-tagging algorithm adopted in our analysis is the one defined as 'loose tag- 
ging' in PGS [69]. Needless to say, a fully rigorous treatment of this issue requires a dedicated experimental 
analysis. 
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(c) Transverse sphericity St > 0.15. 7 

As this process involves several jets, in order to reduce the uncertainty associated with 
jet momentum resolution we use here the m,T2 defined in terms of the total transverse mass, 
rriTv, of the 4 jets in the gluino decay g — > ti\ — > bqqb\\ , see eq. (2.10). Furthermore, in 
the procedure to compute m^y, we do not use the invariant mass of each jet recorded in 
the detector, but simply set rrij = for non-6 jets, and rrij = rrib for 6-tagged jets. This 
means that, for a given trial chargino mass, only the jet transverse momenta are used for 
the calculation of jxiti- Quite often, this approach turns out to be useful for identifying the 
endpoint when the decay product of each mother particle involves more than one (or two) 
jets. Still, the calculation of rriT2 is subject to combinatoric uncertainty due to correctly 
pairing the 4 b jets with the two reconstructed W bosons. To perform this pairing, we use 
the so-called mTGen pairing scheme [70]. 8 

With the events passing the above selection criteria, 9 we have constructed the mT2( m x) 
distribution, with trial chargino mass m x £ [0, 200] GeV. The resulting density plot is 
shown in fig. 8 (upper- left panel). The kink structure can be roughly located by inspecting 
the uppermost dark lines (indicating the largest density of events) , if one neglects the more 
sparse (lighter in color) contributions from background events. 

A quantitative determination of the rriT2 maximum line can be carried out as follows. 
First, one fits the endpoint of the rriT2 distributions obtained at fixed m x . Examples of 
these distributions, and their endpoint fits to a two-segments line, are reported in fig. 8 
(lower panels) for the values m x = {0, 100, 197} GeV. Here and henceforth (figs. 8-11), the 
two-segments line is parameterized as 

y = pi{x-p ) +P2X + P3 , fora:<po, ^^ 

y = p 2 x + p 3 , for x > po . 

The resulting endpoint values of rriT2 for the various values of m x , namely m T 1 2 X ( m x)' are 
shown in fig. 8 (upper-right panel), with bars representing the statistical error. To identify 
the kink position, we then fit those endpoint values with appropriate fitting functions. 



7 We adopt the same definition as [44], namely 

* ■ >^h • ("' 

where Ai,2 are the eigenvalues of the 2x2 sphericity tensor Sij — ^kPip'j, with k running on all the 
visible particles, and i,j = {1,2} labeling the momentum transverse component. SUSY events tend to be 
spherical (St ~ 1) as opposed to back-to-back. 

8 In the present example, the rriTGcn scheme works as follows. Assuming that only one of the (three) di-jet 
pairings among q jets fulfills the Mw requirement, there are still six possible ways of forming two b di-jets 
and of assigning them to the two reconstructed W bosons. For each of these combinations, we calculate the 
mT2(m x ) value for fixed m x , and pick the combination corresponding to the minimum mT2(m x ). As m x 
value we took 100 GeV and we checked that the pairing choice does not change when reasonably varying 
m x . 

9 In particular, we checked by a parton-level analysis that the adopted wiTGen pairing scheme is effective 
in reducing the jet combinatoric error, while not significantly affecting the final txiti distribution with 
respect to the true pairing. 
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Figure 8: Upper-left panel: rriTiijn^) density plot for the events of step 1, LS scenario. The color 
code on the right of the plot represents the number of events described by each line; Upper-right 
panel: fit to the corresponding maximum TOT 2 (m x ); Lower panels: m,T2 distribution for the events 
of step 1, LS scenario, with m x = {0, 100, 197} GeV. 



To this end, one can use the known analytic form of mT2(m x ) [^Ij 54] for events with 
m,T2(m x = m~±) = rrig and negligible upstream transverse momentum: 



m T 2(m x ,m v ) 



m\ — m 2 _+ + m?, 
9 Xi v _ 

2rrig 



+ 



ml + m?+ 

9 Xi 



mi 



+ 4m? 



m: 



2 

m + 

Xi 



2m g 



,(4.3) 



where my denotes the total invariant (or transverse) mass of the 2b + 2q jets in the gluino 
decay. From this m,T2 functional form, one can easily notice that m™2 X = mT2{ rn v = m™ ax ) 
for m x > m-±, while rrij^ = mT2(my = rriy m ) for m x < m~±, hence m™^ shows a kink 
at m x = m~±. Generically, m™ ax < m^ — m~± and m.™ 111 > 0, and one easily finds 
m -± + m y for m?? ax = mg — m-± . 

Xl A V U X -^ 

The rriT2 endpoint values in fig. 8 for m x > m~± indicate that they can be well 
described by a straight line, implying that my ax is close to m~ g — m~± . Inspired by the 



T2 
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analytic form (4.3), we then use the following fitting functions to find the kink position: 



max 



= Mi + A /M? + ml if m y < m~± 



,u T2 - ±v±i_ t y ivj. 2 T m- x 11 in x ^ m-± , ^ 

m™| x = M3 + m x if ?n x > m~± , 



where Mi (i = 1,2,3) are fitting parameters. The best-fit curves are shown in fig. 8 
(upper-right panel), together with the corresponding values of Mj. The resulting gluino 
and chargino masses are given by rn~± = 109(17) GeV and rrig = 395(16) GeV, which are 
fully compatible with the true values in table 1. We also find that the m^^-curves in fig. 
8, which describe data at detector level, are reasonably close to the curves obtained at 
parton level. 



m-o 



Step 2: For t\t\ -^2b+4Q+i ) T events, determine the mass differences m^ 

and m~± — m y o from the endpoints of the transverse-mass distributions rriT,bqq an d iriT,qq, 

respectively 

To select these events, we need to require 2 jets tagged as b, 4 q jets, and no leptons. 
The channel of t\ pair production is subject to SM backgrounds coming mostly from tt or 
WW + 26 production. These backgrounds can be eliminated efficiently by the following 
selection cuts 

(a) pti,2 > 50, 25 GeV on the 2 b jets; 

(b) pri.2,3,4 > 50, 25, 20, 10 GeV on the 4 q jets; 

(c) $r> 100 GeV; 

(d) Transverse sphericity St > 0.15. 

For the events passing the above cuts (a)-(d), we calculated, event by event, all possible 
di-jet transverse mass values out of the 4 q jets and picked the minimum value among them. 
Here we use the di-jet transverse mass thtjj, rather than the di-jet invariant mass mjj, 
as the endpoint of m,jj is severely smeared. The corresponding distribution is shown in 
fig. 9 (upper-left panel). Fitting the endpoint region to a two-segments line, we obtain 
mS^j = m-± — m^o = 51.8(1.7) GeV, corresponding to m^o ~ 58 GeV, a good estimate 
of the true value in table 1. A completely analogous procedure can be adopted for the 
transverse mass of the bqq' system, and the resulting distribution is shown in fig. 9 (upper- 
right panel). The fitted endpoint m^ of this distribution allows to determine the t\ mass 
through the following formula 



m h 



( m max )2 + m 2 _ m 2 

K — T 2 7% — - . ( 4 - 5 ) 



implying m^ ~ 206 GeV, again consistent with the true value in table 1. 



in 



10 Note that, if the relation m? ± = m t - m^o is fulfilled (as is our case within about (25 GeV) 2 ), eq. (4.5) 
simplifies to m™ ^ = m^ — m^o . 
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Figure 9: Upper-left panel: transverse mass distribution m,T,qq (step 2, LS scenario); Upper-right 
panel: bqq-jct transverse mass distribution (step 2, LS scenario); Lower panel: invariant mass 
distribution of the £ + £~ pair (step 3, LS scenario). 



Step 3: For \\ X2 ~ * ^^ + ^ + i>T events, determine the mass difference m^o — m^o 
from the endpoint of the invariant-mass distribution nun 

To select these events, we require 2 leptons (either e or fj,) of same flavor and opposite 
charge, one additional lepton of different flavor (with veto on hadronically decaying taus), 
and no jets. To clean up the event selection from possible backgrounds, we also impose the 
following cuts: 

(a) pri,2 > 20, 15 GeV on the 2 leptons of same flavor; 

(b) pt > 15 GeV on the lepton of different flavor; 

(c) $ T > 50 GeV; 

(d) Transverse sphericity St > 0.15. 

The invariant mass distribution constructed with the 2 leptons of same flavor is shown 
in fig. 9 (lower panel). The fact that, in general, the lepton signal is quite accurately 
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reconstructed explains the clear end point structure of the invariant mass distribution at 
50 GeV (even with a limited number of events). This endpoint value is in agreement with 



the expected m^o 



A 2 






51 GeV from table 1. 



4.2 Analysis: HS scenario 

Step 1: Construct the niT2 variable for gg — )• 4 b (+ 2 ^) + $t decay events, thereby 
measuring m,g and m^o 




o 540 

U 

lp 520 
E 

500 




7 


X 2 /ndf 0.5501/9 
M, 282 ±16.33 
M, 120.3± 18.53 
M 3 330.1 ±2.059 




480 
460 
440 
420 
400* 



















100 120 140 160 180 200 

n\ [GeV] 



I4.K6/7 
S 10.6 ±4.6 
1. SIR* 0.186 
0.221+0.054 
138.9+ 11.0 




100 200 300 400 500 



Figure 10: Upper-left panel: m T2 (m x ) density plot for the events of step 1, HS scenario. The color 
code on the right of the plot represents the number of events described by each line; Upper-right 
panel: fit to the corresponding maximum rriT2(m x ); Lower panels: m,T2 distribution for the events 
of step 1, HS scenario, with m x = {0, 100, 180} GeV. 



Within gg production, the decay g — > ^bb, followed by x® ~^ X?7 is the dominant channel. 
Since this channel is the only possible one giving 4 jets, 2 7 and no lepton, it may not be 
necessary to require b tagging. Indeed, we select the above event content without b tagging. 
We however impose the following cuts as well: 



(a) pn.2,3,4 > 120 > 70 ' 30 > 20 GeV on the 4 J ets ; 

(b) pti,2 > 50, 20 GeV on the hard photons; 
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(c) $ T > 100 GeV; 

(d) Transverse sphericity St > 0.15. 

On the selected events, we calculate mr2(%) with the 4 jets according to the rriTGen 
pairing scheme, and treating the 2 7 as part of the $t- Here we use the di-jet invariant mass 
rrijj, rather than the transverse mass, to compute mn, as the resulting rriT2 distribution 
shows a reasonably clean endpoint structure. 

The rriT2{fnx) density plot, displayed in fig. 10 (upper- left panel), shows a substantial 
number of spurious events contributing to the maximum mT2 line region, impairing a 'by 
eye' determination of the kink. Still, we can quantify the kink position by following the same 
strategy as step 1, LS scenario, namely fitting first the endpoint of the rriT2 distributions 
obtained at fixed m x values. Fig. 10 (lower panels) shows those distributions and fitted 
endpoints for m x = {0, 100, 180}. Fitting those endpoints (obtained for various values of 
m x ) to the function in eq. (4.4) (see fig. 10, upper-right panel), we get the kink position at 
m x o = 126(16) GeV and m g = 456(15) GeV, which are in agreement with the true values 
in table 1. (The best-fit values of M, are also reported in the figure.) Comparing to the 
parton level result, we find the kink to be not as sharp as expected. The main reason for 
this is that a large portion of the events with small rrijj, that provide the true to™| x for 
m x < m x o, are eliminated by the selection cuts. 

Step 2: For X\ X2 ~^ f- + ^~ + & + i>T events, determine the mass difference m x o — m x o 
from the endpoint of the invariant-mass distribution mu 

This step is completely analogous to step 3 of the LS scenario. We require two leptons (e 
or /i) of same flavor and opposite charge, one lepton of different flavor n , and no jets. To 
remove possible backgrounds, we also impose the following cuts: 

(a) pti,2 > 20, 10 GeV on the 2 leptons of same flavor; 

(b) pt > 10 GeV on the lepton of different flavor; 

(c) $ T > 50 GeV. 

The invariant mass distribution constructed with the 2 leptons of same flavor is shown in 
fig. 11 (upper-left panel). Again, the accuracy with which lepton signals are in general 
reconstructed explains the clear endpoint structure of the invariant mass distribution at 60 
GeV. This endpoint agrees with the expected value m^o — m^a = 59 GeV from table 1. 

A 2 Al 

Step 3: For X\ X2 ~ * ®Q + ^ + i*T events, determine the mass difference tti~± — m x o 
from the endpoint of the invariant-mass distribution m qq (and again m x o — m x o from mu) 

We select these events by requiring 2 q jets and 2 leptons (e or //) of same flavor and 
opposite charge. We also impose the following cuts : 



At variance with step 3 of the LS scenario, here we also allow this one lepton to be a hadronically 
decaying tau in order to have larger statistics. 
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Figure 11: Upper-left panel: invariant mass distribution of the £ + £~ pair (step 2, HS scenario); 
Upper-right panel: di-jet invariant mass distribution from Xi (step 3, HS scenario); Lower panel: 
invariant mass distribution of the £ + £~ pair (step 3, HS scenario). 



(a) p T1>2 > 20, 10 GeV on the 2 jets; 

(b) pti >2 > 20, 10 GeV on the 2 leptons; 

(c) $ T > 50 GeV. 

The di-jet invariant mass distribution is reported in fig. 11 (upper-right panel). As can be 
seen by the length of its tail, the background contribution is quite severe for this channel 

even after the above cuts. The fitted endpoint allows to estimate m~± — m^o to be 78(11) 

Xi X-l v ' 

GeV, whose central value is 35 % larger than the true value 58 GeV. So, because of this 
systematics, the accuracy in the determination of rn-± turns out to be not as good as for 
the other mass determinations in the HS scenario. 

In addition, we can check m^o — m^o as obtained in step 2 by constructing the invari- 
ant mass of the two leptons. Fig. 11 (lower panel) displays the corresponding distribution, 
showing again the quite clear endpoint at 60 GeV, though it is more smeared by back- 
grounds than the distribution in step 2. 
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5. Conclusions and Outlook 

We have explored, in the concrete example of Yukawa-unified SUSY GUTs, to which extent 
a quantitative determination of the lightest part of the SUSY spectrum is practicable with 
LHC data. Specifically, we have considered two representative scenarios, both viable in 
the light of all existing data, but characterized by some notable differences in the SUSY 
spectra (summarized in table 1), that arguably can be told apart only via their direct 
measurement. 

We have elaborated a mass-determination strategy based on the observation of kinks 
in the kinematic variable rriT2i namely the so-called TO*r2-kink method. This method is 
especially suitable for our purposes, as it does not rely on the presence of long decay 
chains, which are not achievable at the LHC within our considered classes of models. 

We have studied in detail our strategy on 100 
fb _1 of data collected at 14 TeV p-on-p collisions, 
using Pythia 6. We have also simulated detector- 
level effects via PGS4, with the aim of conveying a 
hopefully realistic idea of how large a degradation 
is expected in real signals. 

We were able to determine the masses of the 
gluino, of the lightest chargino, of the first two neu- 
tralinos for both scenarios considered, and also the 
mass of the lightest stop for the scenario where this 
mass is below the gluino's. Our results are summa- 
rized in table 5, as obtained through a global fit to 
the outcomes (masses or functions thereof) of the 
various steps described in section 4 for either sce- 
nario. The results in table 5 display a good agree- 
ment with the corresponding true values in table 1, 
a statistical error always around 20 GeV and a systematics somewhat larger for the HS 
scenario. 

Our worked example of Yukawa-unified SUSY GUTs may offer a useful playground for 
dealing with other theories which predict similar patterns of SUSY spectra. 

Of course, a number of issues are left open by the present study. First, although we 
have been focusing on mass determinations, it is also of crucial importance to determine 
the spin of any SUSY particle produced at the LHC. To this aim, we just note here that the 
approach of ref. [71] , whereby the invisible particle momenta are reconstructed with m,T2 in 
order to determine the mother-particle spin, may be applied to the decay of pair-produced 
gluinos in the scenarios of [10, 16], from which the gluino spin might be read off. 

Further issues are more specifically concerned with mass measurements and the preci- 
sion achieved in their determination. For example, a first relevant question is that of the 
luminosity needed for the discovery of either channel. A rough extrapolation of our results 
indicates that about 5 to 10 fb _1 of data would be sufficient to obtain a statistical error 
on the mass determinations of about 50 GeV. This amount of data would therefore suffice 
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66 ±16 
126 ±16 



Table 5: Mass determinations within 
our strategy, to be compared with table 
1. 
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for a 5cr discovery of at least the heavier among the particles listed in Table 5 for either 
scenario. 

We note as well that we did not try and exploit here the approach of ref. [72], that 
has been proposed as a strategy to make the endpoint structure more prominent, and 
that might be able to reduce the systematic errors in the identification of the endpoint 
position of the kinematic variables considered in our analysis. Pursuing these, more refined 
strategies could also offer a handle to address the case of lower needed luminosities, where 
the signal reach may be increased by also using more optimized cuts, or allow to implement 
our strategy to data from Tevatron Run II and from the LHC running at center-of-mass 
energies below the design one. 
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